Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1252 1758 (N55W27 L=163)  1200    220        1 (1150/1723-1347/1794)  10.1 W14-W42 N52-N58
      2  404 1462 (N32E50 L= 86)  1718    334        2 ( 324/1351- 518/1576) -23.0 E69-E36 N24-N39
      3 1164 1548 (N38W12 L=148)   578    140        0 (1111/1507-1236/1575) -32.0 W07-W16 N34-N39
      4 1343 1505 (N34W25 L=161)   336     55        0 (1308/1488-1379/1520)   9.5 W22-W28 N33-N35
      5  865 1249 (N18E10 L=126)   604    203        0 ( 820/1171- 929/1304)  58.0 E12-E05 N13-N21
      6  966 1267 (N19E03 L=133)   208     60        1 ( 944/1242- 989/1293) -33.7 E04-E01 N17-N20
      7  155 1222 (N15E74 L= 62)  1730    107        4 ( 126/1180- 186/1255)  39.8 E83-E67 N11-N16
      8  995 1226 (N16E02 L=135)   215     33        1 ( 969/1208-1020/1243)   0.0 E03-W00 N15-N17
      9  544  979 (N01E31 L=105)  1390    412        0 ( 463/ 836- 663/1096)  47.0 E37-E22 S07-N08
     10  455  661 (S19E41 L= 95)  1298    235        1 ( 421/ 557- 510/ 737)  53.1 E47-E35 S26-S14
     11 1201  490 (S31W12 L=148)  4079    702       10 ( 943/ 330-1434/ 700)  37.4 E07-W27 S44-S17
     12  844  513 (S29E14 L=123)  1710    543        1 ( 680/ 381-1050/ 691)  35.0 E29-W01 S39-S17
  

    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