Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  559 1699 (N21E54 L= 39)  2054    288        3 ( 505/1609- 623/1807) -26.6 E69-E43 N16-N29
      2 1278 1272 (N13W08 L=100)   435    171        0 (1244/1208-1325/1346) -32.0 W03-W11 N10-N16
      3 1385 1268 (N16W14 L=107)   345     74        0 (1352/1245-1427/1286)  11.3 W11-W16 N14-N17
      4 1820  731 (N01W63 L=156)   426     42        1 (1800/ 707-1834/ 761) -26.6 W60-W65 N00-N02
      5 1279  487 (S29W34 L=126)  1825    340        1 (1161/ 417-1422/ 601)  69.8 W29-W37 S37-S18
      6 1582  481 (S19W56 L=148)   316     66        0 (1547/ 469-1616/ 490)  26.6 W52-W58 S20-S17
      7 1532  444 (S23W54 L=147)   635     81        2 (1498/ 408-1565/ 476) -45.0 W50-W57 S25-S21
      8 1036  359 (S45W26 L=119)  2331    281        9 ( 910/ 296-1146/ 430)   2.4 W13-W37 S48-S41
  

    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