Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1650 1276 (N28W34 L= 33)   314     49        1 (1633/1251-1671/1305) -33.7 W32-W35 N26-N29
      2  147 1166 (S15E69 L=290)   638    103        3 ( 129/1136- 164/1205) -78.7 E71-E66 S17-S12
      3 1752 1105 (N20W45 L= 44)   841    219        1 (1688/1076-1825/1131)  13.0 W39-W52 N18-N22
      4 1282 1055 (N04W14 L= 13)   198     41        0 (1263/1027-1303/1084)  90.0 W12-W14 N02-N06
      5 1595  955 (N07W36 L= 35)   330    120        0 (1561/ 904-1621/1020)  82.9 W34-W37 N03-N11
      6  392  943 (S24E39 L=321)   305     66        1 ( 371/ 907- 412/ 979)  78.7 E40-E37 S25-S20
      7 1281  666 (S17W25 L= 24)   206    119        0 (1216/ 652-1332/ 684)  29.7 W20-W27 S18-S14
      8 1533  586 (S15W45 L= 44)   488     92        2 (1491/ 555-1565/ 620)  53.1 W41-W47 S16-S12
      9  912  457 (S41W10 L=  9)   829    154        2 ( 863/ 400- 973/ 523) -22.6 W04-W16 S43-S36
     10 1071  408 (S39W23 L= 22)   748    147        1 ( 999/ 382-1143/ 441)  15.3 W18-W29 S41-S37
     11 1562  398 (S24W58 L= 57)  1135    267        1 (1487/ 379-1637/ 423)  15.9 W50-W65 S26-S21
  

    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