Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1341 1735 (N51W39 L=281)   631     80        0 (1305/1699-1362/1760) -90.0 W36-W42 N47-N53
      2 1848 1327 (N16W71 L=313)  1147    121        0 (1792/1289-1898/1367) -20.6 W64-W81 N13-N19
      3  492 1155 (N14E35 L=207)   319     64        2 ( 458/1130- 515/1183) -14.0 E37-E33 N12-N16
      4  332  935 (N01E48 L=193)   301     52        1 ( 317/ 904- 351/ 963) -56.3 E49-E46 S01-N02
      5  710  871 (S04E20 L=221)   669    271        1 ( 592/ 827- 824/ 933) -18.4 E28-E13 S06-N00
      6 1116  766 (S12W04 L=246)   202     77        0 (1085/ 751-1161/ 777)  11.3 W02-W07 S12-S11
      7  930  697 (S16E08 L=234)   238     61        0 ( 907/ 668- 957/ 725)  53.1 E09-E05 S17-S13
      8 1577  629 (S24W38 L=280)   999    199        3 (1507/ 554-1641/ 703)  45.0 W33-W42 S29-S19
      9  310  631 (S19E57 L=184)  1129    133        2 ( 283/ 573- 353/ 684)  45.0 E64-E51 S23-S15
     10  733  531 (S26E23 L=219)   365    153        0 ( 690/ 478- 796/ 571)  33.7 E26-E17 S29-S23
  

    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