Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  674 1380 (N35E17 L=168)   235     55        0 ( 646/1355- 694/1397)   0.0 E18-E14 N33-N36
      2  533 1362 (N37E28 L=157)   251     53        1 ( 514/1332- 554/1390)  53.1 E29-E26 N34-N38
      3  197 1189 (N30E59 L=126)  1673    249        2 ( 156/1108- 242/1257) -51.3 E64-E54 N24-N34
      4  417  858 (N07E42 L=143)  1587    542        3 ( 298/ 748- 500/1008) -62.1 E50-E37 N00-N17
      5 1254  931 (S03W12 L=197)   418     86        2 (1220/ 903-1291/ 975) -68.2 W09-W13 S05-N00
      6 1323  921 (S05W16 L=201)   235     40        2 (1302/ 899-1353/ 944) -45.0 W14-W17 S06-S03
      7  446  487 (S16E52 L=133)  1259    205        4 ( 403/ 444- 495/ 539) -56.3 E56-E47 S19-S12
  

    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