Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1141 1576 (N35E05 L=178)  2573    518        9 ( 892/1453-1405/1714)  -7.9 E27-W16 N30-N39
      2  503 1556 (N20E48 L=135)   746    313        2 ( 472/1449- 535/1651) -42.7 E55-E42 N13-N25
      3 1463 1436 (N33W21 L=204)   244     51        1 (1441/1417-1489/1453)  33.7 W19-W23 N31-N34
      4  538 1191 (S01E33 L=150)   408     79        1 ( 495/1159- 568/1228)  45.0 E35-E31 S03-N01
      5  980  882 (S11W00 L=183)   186     46        0 ( 959/ 856-1004/ 908) -33.7 E01-W02 S12-S09
      6 1152  850 (S09W11 L=194)   741    204        4 (1057/ 817-1257/ 908)  32.5 W05-W16 S12-S04
      7 1506  566 (S19W42 L=225)   310     66        1 (1482/ 549-1534/ 585)  56.3 W39-W44 S19-S16
      8 1137  513 (S30W19 L=202)   465    102        2 (1095/ 488-1183/ 548)  -8.1 W15-W22 S31-S28
  

    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