Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  557 1340 (N35E26 L=169)   253     57        1 ( 537/1312- 581/1362)  76.0 E27-E24 N32-N36
      2  167 1105 (N24E63 L=133)   504    132        1 ( 140/1070- 190/1139) -45.0 E66-E60 N22-N26
      3 1069  889 (S02E00 L=196)   244     51        2 (1043/ 859-1089/ 920) -90.0 E01-W01 S03-N00
      4 1160  869 (S05W05 L=201)   395     91        2 (1111/ 830-1205/ 911) -56.3 W02-W07 S07-S01
      5  313  806 (N05E53 L=143)  1119    183        1 ( 253/ 741- 369/ 889) -59.0 E57-E48 N00-N10
      6  997  801 (S06E06 L=190)   239     53        0 ( 964/ 775-1028/ 820)  11.3 E08-E03 S06-S04
      7 1152  770 (S10W03 L=198)   372    110        0 (1107/ 736-1194/ 794) -36.9 W00-W04 S13-S08
      8 1062  685 (S14E04 L=191)   247     57        3 (1041/ 650-1080/ 728)  53.1 E06-E03 S15-S11
      9  349  461 (S18E64 L=132)  1878    314        5 ( 299/ 425- 420/ 510)   0.0 E70-E55 S20-S13
  

    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