Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1186 1751 (N55W19 L=167)   724    163        2 (1121/1724-1269/1776)   6.3 W11-W29 N52-N56
      2  968 1557 (N39E04 L=145)   638    164        2 ( 888/1523-1043/1601)  -9.5 E10-W02 N36-N42
      3  324 1450 (N30E59 L= 89)  2024    270        2 ( 277/1346- 382/1576) -21.3 E85-E48 N23-N37
      4 1184 1498 (N34W12 L=161)   251     50        0 (1155/1479-1215/1515)  21.8 W10-W15 N32-N35
      5  704 1257 (N18E21 L=128)   704    229        1 ( 638/1183- 774/1308)  29.1 E25-E16 N13-N21
      6  363  971 (S00E45 L=103)  1655    338        3 ( 299/ 849- 465/1105)  49.1 E51-E37 S07-N08
      7 1136  878 (S05W07 L=155)   327    143        0 (1078/ 837-1192/ 910) -29.7 W03-W10 S07-S03
      8  316  670 (S19E54 L= 94)  1617    203        3 ( 290/ 582- 365/ 736)  47.7 E60-E48 S25-S14
      9 1608  693 (S18W41 L=189)   410    154        0 (1533/ 676-1667/ 724)  10.3 W34-W46 S19-S16
     10  737  572 (S25E20 L=128)  1480    457        2 ( 583/ 429- 904/ 720)  34.6 E36-E07 S36-S15
     11 1021  481 (S32E01 L=148)  3469    671       10 ( 778/ 344-1248/ 692)  36.3 E21-W14 S42-S17
     12 1624  645 (S22W43 L=192)   311     61        2 (1594/ 624-1654/ 668)  11.3 W40-W46 S23-S20
  

    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