Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  467 1588 (N16E52 L=218)   475     73        0 ( 441/1565- 492/1617) -15.9 E55-E48 N15-N18
      2 1328 1533 (N30W08 L=277)   473    156        1 (1292/1465-1360/1605)  90.0 W06-W09 N25-N36
      3 1495 1465 (N30W21 L=290)  1246    370        0 (1389/1314-1595/1597) -35.8 W11-W29 N22-N37
      4  427 1528 (N12E52 L=218)   357     60        2 ( 407/1505- 454/1555)  33.7 E54-E49 N10-N13
      5  592 1532 (N15E39 L=231)   467     91        1 ( 560/1497- 634/1554)  53.1 E41-E36 N12-N16
      6  919 1452 (N16E15 L=254)   693    169        1 ( 858/1423-1002/1500)  -5.7 E20-E10 N14-N17
      7 1567 1093 (N09W31 L=301)  1462    322        1 (1538/ 950-1617/1233) -76.8 W29-W34 N01-N18
      8 1691 1103 (N13W40 L=310)   376    116        0 (1673/1049-1705/1154) -51.3 W38-W42 N10-N15
      9 1825  954 (N09W55 L=325)  1092    180        0 (1750/ 903-1870/1014)  41.6 W48-W60 N04-N12
     10 1620  917 (N01W39 L=308)   294     72        0 (1583/ 891-1654/ 937)   0.0 W35-W41 N00-N01
     11 1376  740 (S15W28 L=297)   490    248        0 (1274/ 694-1485/ 787)  39.3 W22-W33 S19-S10
     12 1576  599 (S18W47 L=316)   723    122        0 (1525/ 563-1624/ 646)  66.8 W43-W48 S20-S13
     13 1725  601 (S13W60 L=329)   536     90        2 (1700/ 559-1751/ 639) -78.7 W57-W61 S15-S10
     14 1057  571 (S33W13 L=283)   226    110        0 ( 996/ 563-1106/ 580)  23.2 W09-W16 S34-S30
     15 1721  521 (S17W65 L=335)  1033    165        2 (1648/ 490-1790/ 564)  39.3 W59-W72 S21-S12
     16 1369  456 (S31W39 L=309)   663    104        0 (1319/ 426-1421/ 492)  66.8 W36-W41 S34-S27
     17 1244  421 (S37W33 L=302)   284     15        1 (1223/ 402-1263/ 443) -33.7 W30-W34 S38-S35
     18 1464  298 (S38W61 L=331)   525     57        0 (1437/ 282-1495/ 320)   0.0 W56-W65 S38-S36
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y