Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1119 1610 (N34E11 L= 76)  1108    200        2 (1077/1527-1154/1691) -42.3 E16-E05 N28-N39
      2  497 1476 (N11E46 L= 41)   482    119        0 ( 479/1423- 519/1540) -66.8 E47-E44 N07-N14
      3  740 1426 (N13E27 L= 60)   315     98        1 ( 699/1393- 781/1456)  -8.1 E30-E23 N11-N14
      4 1420 1237 (N19W18 L=105)   898    217        3 (1365/1169-1472/1334)  72.9 W16-W20 N13-N26
      5 1865 1286 (N36W62 L=149)   829     95        0 (1826/1248-1891/1312)   4.4 W55-W68 N33-N37
      6 1673 1107 (N18W39 L=126)   849    148        1 (1641/1039-1703/1174) -49.4 W35-W42 N14-N21
      7 1086 1128 (N04W01 L= 88)   275     73        2 (1046/1107-1122/1151)  63.4 W00-W02 N02-N06
      8 1803  912 (N10W56 L=143)  1187    246        2 (1761/ 798-1833/ 974) -73.3 W51-W59 N04-N14
      9  456  795 (S30E32 L= 55)   330     76        3 ( 429/ 772- 481/ 830)  51.3 E34-E30 S32-S27
     10  345  782 (S34E43 L= 44)   813    181        3 ( 269/ 743- 419/ 828)  -7.1 E51-E35 S35-S32
     11 1554  629 (S12W44 L=131)   379    134        0 (1507/ 586-1584/ 684)  81.9 W42-W45 S15-S08
     12 1294  538 (S24W31 L=118)   249     85        1 (1270/ 505-1321/ 576) -23.2 W27-W34 S25-S22
  

    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