Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  368 1595 (N11E64 L= 16)  1420    285        4 ( 318/1485- 421/1665) -81.3 E67-E61 N02-N15
      2  801 1607 (N21E31 L= 49)   296    165        0 ( 756/1555- 826/1660) -29.1 E36-E27 N18-N23
      3 1230 1374 (N19W02 L= 81)   807    287        0 (1206/1252-1261/1504) -68.2 E01-W06 N12-N27
      4 1576 1374 (N30W24 L=104)   301     54        1 (1548/1350-1600/1397) -11.3 W21-W26 N28-N31
      5  234 1316 (S09E63 L= 17)   908    173        4 ( 205/1285- 276/1357)  -5.7 E67-E57 S10-S06
      6  743 1281 (N00E23 L= 57)   405    134        0 ( 681/1244- 813/1304)   0.0 E26-E17 S01-N01
      7 1422 1226 (N16W17 L= 97)   463    124        0 (1403/1169-1443/1294) -60.3 W14-W18 N12-N19
      8 1238  929 (S06W14 L= 94)   482     99        1 (1215/ 871-1258/ 984) -74.1 W12-W15 S09-S02
      9  168  961 (S31E59 L= 20)  1562    412        3 (  88/ 938- 246/ 981)  10.0 E68-E51 S32-S28
     10  665  891 (S23E17 L= 62)   694    128        3 ( 614/ 835- 711/ 944)  63.4 E19-E15 S27-S19
     11 1825  865 (N10W57 L=137)   781    151        1 (1796/ 808-1842/ 944) -54.5 W54-W60 N06-N13
     12  420  761 (S38E34 L= 46)  1377    307        3 ( 317/ 726- 540/ 811)  19.3 E43-E23 S42-S33
     13  913  620 (S32W06 L= 85)   957    254        2 ( 831/ 540-1017/ 727)  76.0 W02-W08 S39-S23
     14 1838  683 (N01W68 L=148)   999    269        1 (1780/ 663-1892/ 705)  17.4 W61-W77 S01-N04
     15 1253  610 (S22W26 L=106)   466    143        0 (1177/ 578-1308/ 663)  60.9 W22-W28 S26-S17
     16 1530  597 (S14W44 L=124)   351    117        0 (1493/ 579-1570/ 627)   0.0 W40-W47 S15-S12
     17 1274  467 (S29W34 L=114)   381     58        1 (1238/ 452-1314/ 484)  21.8 W31-W36 S31-S27
     18 1424  318 (S32W53 L=133)   443     54        1 (1393/ 305-1456/ 339)  21.8 W51-W56 S33-S30
  

    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