Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  360 1489 (N17E57 L=283)  1827    329        4 ( 270/1375- 439/1628)  77.5 E63-E53 N08-N26
      2  382 1409 (N12E52 L=288)   577    104        4 ( 347/1354- 421/1460)  66.8 E53-E49 N08-N15
      3 1638 1150 (N17W38 L= 18)   668    264        1 (1599/1016-1667/1263) -66.8 W34-W40 N09-N23
      4  153 1173 (S06E71 L=269)   905    224        1 ( 121/1153- 180/1205)  -8.1 E80-E65 S06-S04
      5 1289 1123 (N10W14 L=354)   624    168        0 (1252/1076-1355/1202) -40.6 W11-W19 N07-N13
      6  747 1157 (N03E19 L=321)   257     32        2 ( 721/1136- 774/1179)  18.4 E20-E17 N01-N04
      7  165 1026 (S14E65 L=275)   895    124        3 ( 148/ 982- 192/1084) -74.1 E67-E61 S17-S10
      8 1608  536 (S20W51 L= 31)   354    128        0 (1580/ 513-1645/ 565) -14.0 W47-W55 S21-S18
      9 1728  540 (S18W63 L= 43)   517    153        0 (1686/ 526-1767/ 558)  18.4 W58-W67 S19-S16
  

    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