Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1374 1623 (N30W43 L=189)  1507    352        4 (1226/1584-1526/1644) -19.2 W31-W54 N24-N35
      2 1690 1478 (N12W58 L=204)   806    132        1 (1644/1437-1728/1529) -74.1 W55-W60 N08-N15
      3 1159 1401 (N24W19 L=164)   750    265        1 (1036/1368-1273/1444) -34.7 W11-W24 N19-N29
      4 1824 1267 (S03W60 L=206)   476     43        1 (1802/1238-1843/1291)  33.7 W57-W63 S04-S01
      5 1848 1072 (S15W56 L=202)  2093    264        5 (1757/ 979-1908/1184)  15.3 W45-W68 S19-S10
      6  638 1085 (N20E21 L=125)  1168    359        2 ( 494/1019- 789/1176) -43.2 E29-E13 N12-N27
      7 1792 1125 (S10W52 L=198)   947    116        2 (1754/1083-1843/1176) -58.0 W49-W56 S13-S05
      8 1919  983 (S23W65 L=211)   708     74        2 (1893/ 945-1942/1023)  14.0 W60-W70 S24-S20
      9 1040  890 (S01E02 L=143)   336    129        0 (1007/ 846-1077/ 926)  21.8 E05-W00 S03-N01
     10  641  862 (N07E26 L=120)   234     86        1 ( 599/ 836- 666/ 880)   0.0 E29-E24 N05-N07
     11 1472  694 (S24W18 L=164)  1840    433        3 (1367/ 541-1579/ 856)  31.8 W07-W28 S31-S16
     12 1776  784 (S29W45 L=190)   734    106        2 (1749/ 735-1799/ 835)  71.6 W41-W47 S32-S26
     13  780  638 (S09E24 L=122)   254    144        0 ( 767/ 588- 820/ 698)  81.9 E25-E22 S12-S05
     14  972  603 (S15E14 L=132)   513    107        1 ( 942/ 549-1013/ 651)  71.6 E16-E11 S18-S12
     15  620  465 (S16E41 L=105)   449    104        2 ( 598/ 420- 650/ 510)  63.4 E43-E38 S18-S12
     16  712  421 (S21E37 L=109)   445     98        1 ( 685/ 372- 730/ 475)  54.5 E39-E33 S24-S17
     17 1329  250 (S54E05 L=141)   554    140        0 (1253/ 221-1382/ 270)  -9.5 E12-W00 S55-S52
     18  850  205 (S40E43 L=103)  1454    153        2 ( 787/ 173- 914/ 249)   0.0 E51-E35 S42-S37
  

    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