Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  636 1376 (N37E16 L=158)   844    201        1 ( 567/1323- 714/1413)   4.4 E22-E09 N35-N38
      2  164 1209 (N37E60 L=114)   664    250        1 ( 127/1171- 206/1251) -36.9 E64-E55 N33-N39
      3 1815  990 (S18W51 L=225)  1814    317        3 (1758/ 884-1889/1076)  25.0 W44-W60 S23-S14
      4  368 1006 (N21E43 L=132)   941    195        1 ( 301/ 977- 442/1060) -35.0 E47-E36 N17-N25
      5 1446  788 (S18W18 L=193)   589    186        0 (1382/ 717-1518/ 853)  16.7 W12-W24 S20-S15
      6  214  775 (N10E64 L=111)  1377    189        6 ( 182/ 719- 251/ 819)  56.3 E68-E59 N06-N12
      7  649  559 (S11E35 L=139)  2216    437        4 ( 591/ 386- 740/ 749)  70.8 E38-E28 S23-N00
  

    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