Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  436 1628 (N19E55 L=276)  1306    276        1 ( 405/1582- 488/1701) -26.6 E62-E48 N16-N24
      2  375 1409 (N04E51 L=280)   478    132        2 ( 354/1376- 398/1445) -33.7 E53-E47 N02-N06
      3 1005 1288 (N08E06 L=324)   594    228        2 ( 941/1230-1097/1323)   6.3 E09-W00 N03-N09
      4  376 1083 (S15E42 L=288)   635    182        1 ( 310/1035- 435/1131)  33.7 E47-E38 S18-S11
      5  926 1041 (S08E06 L=325)   230    104        0 ( 873/1026- 985/1065)   0.0 E09-E02 S08-S07
      6 1080  835 (S17W07 L=338)   216     45        1 (1054/ 810-1103/ 862)  76.0 W06-W08 S18-S14
      7  828  775 (S25E07 L=324)   486    159        0 ( 753/ 734- 903/ 825)  45.0 E11-E02 S29-S20
      8  861  203 (S67W18 L=349)   700     92        3 ( 836/ 171- 886/ 235) -15.9 W11-W25 S69-S63
  

    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