Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1029 1724 (N36E18 L= 89)   588     94        2 ( 973/1690-1080/1750)   0.0 E21-E13 N34-N37
      2  383 1547 (N12E56 L= 51)   707     97        4 ( 351/1500- 424/1586)  80.5 E59-E53 N09-N15
      3 1214 1505 (N25W01 L=108)   329     94        2 (1182/1469-1252/1546)  71.6 W00-W02 N22-N28
      4 1517 1512 (N34W22 L=128)  1228    372        0 (1342/1447-1640/1545)   9.5 W08-W32 N29-N37
      5 1053 1262 (N07E03 L=104)   221     55        0 (1022/1246-1083/1281)  45.0 E04-E01 N05-N08
      6  191 1217 (S09E63 L= 44)   706    111        1 ( 162/1162- 222/1266)  56.3 E65-E60 S12-S06
      7  290 1145 (S13E51 L= 56)   332     70        1 ( 252/1127- 322/1165)   8.1 E54-E47 S13-S11
      8 1820 1097 (N16W51 L=158)   763    108        0 (1777/1047-1854/1138) -18.4 W46-W55 N14-N17
      9 1424 1049 (N03W22 L=129)   329     90        0 (1412/ 994-1444/1095) -59.0 W21-W24 N00-N05
     10  879  979 (S13E07 L= 99)   926    525        1 ( 678/ 899-1079/1058)  31.8 E18-W03 S20-S05
     11  296  961 (S23E48 L= 58)   427     56        1 ( 260/ 940- 325/ 985)  23.2 E52-E45 S24-S21
     12 1833  827 (N02W59 L=166)  1568    312        2 (1774/ 679-1868/ 967) -76.8 W56-W62 S07-N10
     13 1120  844 (S15W09 L=116)   176     48        0 (1100/ 817-1142/ 868)  76.0 W08-W10 S17-S13
     14 1283  733 (S18W22 L=129)   244     55        1 (1261/ 705-1307/ 755)  71.6 W20-W23 S19-S16
     15 1117  691 (S24W13 L=120)   247    115        1 (1061/ 672-1166/ 726)   0.0 W08-W16 S24-S22
     16 1789  605 (S11W66 L=172)  1170    183        7 (1730/ 552-1844/ 650)  10.0 W58-W75 S12-S08
     17 1155  607 (S28W18 L=125)   205     48        0 (1129/ 589-1184/ 626)  56.3 W16-W19 S29-S26
  

    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