Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1267 1765 (N52E07 L= 81)   383     78        0 (1233/1748-1305/1782)   0.0 E10-E02 N50-N54
      2 1107 1608 (N34E11 L= 77)  1131    177        3 (1067/1524-1145/1686) -39.3 E16-E05 N28-N39
      3  488 1483 (N11E47 L= 41)   450    121        0 ( 466/1425- 509/1540) -76.0 E48-E45 N06-N14
      4  731 1427 (N13E28 L= 60)   286     92        1 ( 692/1397- 767/1455)  -8.1 E31-E24 N11-N14
      5 1737 1435 (N42W45 L=133)   491     58        1 (1706/1414-1765/1453)  20.6 W41-W49 N40-N44
      6 1010 1325 (N13E08 L= 80)   410    125        5 ( 966/1275-1048/1371) -80.5 E10-E06 N10-N16
      7 1404 1246 (N19W17 L=105)   932    256        2 (1354/1170-1480/1351)  70.3 W15-W20 N13-N27
      8 1857 1288 (N36W61 L=149)  1031    127        0 (1819/1246-1884/1319)  36.9 W54-W66 N32-N38
      9 1662 1106 (N18W38 L=126)  1041    149        2 (1628/1030-1703/1178) -36.9 W34-W42 N13-N21
     10 1067 1130 (N04E00 L= 88)   263     85        2 (1024/1110-1106/1154)  45.0 E02-W01 N01-N05
     11 1797  918 (N11W55 L=143)  1719    264        7 (1749/ 793-1839/1000)  90.0 W50-W59 N03-N16
     12  448  812 (S29E33 L= 55)   428    111        0 ( 416/ 780- 479/ 868)  60.3 E35-E31 S32-S25
     13  332  787 (S34E44 L= 44)  1017    181        5 ( 257/ 745- 405/ 838)  -7.1 E53-E36 S36-S31
     14 1437  622 (S16W36 L=124)   291    108        1 (1397/ 589-1493/ 645)  38.7 W34-W39 S17-S12
     15 1004  598 (S29W10 L= 98)   391    114        1 ( 973/ 549-1048/ 643) -29.7 W07-W14 S30-S26
     16 1273  540 (S25W29 L=117)   419    137        0 (1221/ 509-1330/ 563)   0.0 W25-W34 S25-S23
     17 1616  446 (S20W61 L=149)  1266    245        4 (1535/ 354-1677/ 537) -82.4 W58-W65 S28-S13
     18 1502  422 (S25W53 L=141)   626    165        0 (1470/ 376-1548/ 479) -81.9 W50-W57 S27-S20
     19  932  253 (S56W24 L=112)   621    107        1 ( 878/ 233- 983/ 276)  11.3 W17-W28 S57-S53
  

    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