Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  325 1493 (N07E60 L=137)   526    121        2 ( 307/1444- 346/1533) -56.3 E62-E56 N03-N09
      2 1522 1371 (N25W22 L=220)   372    111        1 (1491/1344-1561/1390)  11.3 W19-W25 N22-N26
      3 1140 1189 (N05W03 L=200)   239     76        0 (1099/1166-1181/1218)  51.3 W00-W04 N02-N07
      4  657 1114 (S10E23 L=174)   361    109        1 ( 635/1051- 673/1170) -63.4 E25-E21 S13-S07
      5  519 1010 (S19E30 L=167)   229     88        0 ( 480/ 981- 547/1045)  59.0 E32-E28 S20-S15
      6 1871  854 (N07W63 L=260)   448     91        1 (1857/ 812-1880/ 905) -59.0 W60-W65 N04-N09
      7 1346  828 (S10W23 L=220)   326    143        2 (1267/ 805-1398/ 852)  20.6 W18-W26 S12-S08
      8 1274  572 (S26W28 L=225)   242     63        0 (1250/ 555-1312/ 593)  26.6 W25-W29 S27-S24
      9  862  215 (S65W22 L=219)   485     71        1 ( 832/ 198- 887/ 241)  90.0 W18-W25 S67-S61
  

    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