Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  998 1451 (N25W03 L= 51)   687    163        1 ( 903/1417-1066/1476) -21.8 E03-W07 N22-N27
      2  316 1327 (N25E48 L=  0)  2175    346        1 ( 191/1233- 444/1441) -25.3 E71-E35 N17-N34
      3  870 1065 (N03E09 L= 39)   184     54        0 ( 840/1052- 899/1077)   0.0 E10-E06 N02-N03
      4 1865  955 (S13W61 L=109)   609    106        1 (1840/ 904-1892/1005)  51.3 W57-W64 S15-S10
      5 1409  764 (S20W22 L= 70)   334     70        0 (1369/ 749-1444/ 790)   0.0 W19-W24 S21-S19
      6 1288  743 (S20W14 L= 62)   191     61        0 (1273/ 714-1309/ 774)  63.4 W12-W15 S22-S18
      7 1815  664 (S31W61 L=109)  1574    172        0 (1766/ 631-1882/ 713)   0.0 W53-W73 S32-S27
      8  474  494 (S27E48 L=  1)  1532    241        0 ( 423/ 384- 528/ 586)  31.0 E59-E39 S34-S21
      9  701  376 (S38E35 L= 14)   288     62        0 ( 668/ 361- 733/ 394)   9.5 E37-E31 S39-S37
  

    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