Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1625 1585 (N18W58 L= 53)  1883    355        5 (1507/1529-1749/1625) -32.3 W49-W68 N10-N22
      2 1735 1256 (S02W51 L= 46)   795    142        1 (1710/1199-1768/1330)  74.1 W48-W54 S06-N01
      3  731 1029 (N14E16 L=339)   250     86        1 ( 692/1001- 780/1047)  -9.5 E19-E13 N12-N14
      4  231  881 (N16E57 L=298)   428     62        1 ( 214/ 849- 250/ 913)  36.9 E58-E54 N14-N17
      5  261  768 (N09E58 L=297)  1515    426        2 ( 149/ 711- 377/ 807) -20.0 E71-E49 N03-N12
      6  948  515 (S21E19 L=336)   350    111        1 ( 917/ 472- 973/ 572)  45.0 E22-E16 S23-S17
      7  552  316 (S24E55 L=300)  1824    362        3 ( 510/ 210- 595/ 407)  55.5 E61-E49 S33-S17
      8 1076  266 (S43E23 L=332)   407    134        1 (1051/ 225-1108/ 316)  31.0 E28-E18 S45-S39
  

    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