Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1328 1771 (N38W49 L= 95)   712    113        0 (1290/1746-1368/1791)   0.0 W44-W52 N36-N40
      2 1411 1770 (N35W54 L=101)   478     48        0 (1388/1743-1431/1790)  11.3 W51-W56 N34-N36
      3 1695 1542 (N14W62 L=108)   579     84        1 (1666/1511-1725/1572) -78.7 W59-W64 N11-N16
      4 1144 1502 (N28W21 L= 68)   431    128        1 (1087/1473-1219/1520) -12.5 W17-W26 N26-N29
      5 1010 1276 (N19W06 L= 53)   342     74        1 ( 965/1252-1052/1295) -31.0 W03-W08 N17-N20
      6  202 1274 (N39E55 L=352)   747    214        0 ( 152/1258- 253/1294) -15.3 E60-E49 N37-N41
      7 1869 1068 (S17W58 L=105)   737    156        4 (1841/1036-1896/1098)  18.4 W55-W62 S18-S14
      8 1778  942 (S21W46 L= 93)   426     84        0 (1751/ 911-1801/ 983) -78.7 W44-W47 S22-S17
      9 1447  687 (S26W16 L= 63)  1102    285        2 (1407/ 557-1512/ 818)  56.3 W12-W23 S34-S19
     10 1141  623 (S21E04 L= 43)   744    235        0 (1099/ 523-1172/ 744)  56.3 E09-W00 S26-S14
     11  850  485 (S22E26 L= 21)   618    146        1 ( 820/ 414- 900/ 543)  84.3 E26-E23 S27-S17
     12  574  298 (S27E54 L=352)   600    131        1 ( 551/ 264- 596/ 337)  32.0 E57-E49 S29-S24
  

    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