Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  447 1709 (N20E64 L=131)  1477    268        2 ( 423/1665- 491/1755) -37.9 E70-E60 N16-N23
      2  667 1456 (N09E33 L=163)   227     22        1 ( 648/1437- 685/1477)  90.0 E34-E31 N07-N10
      3  692 1009 (S16E19 L=177)   561    126        1 ( 628/ 973- 749/1048)  50.2 E22-E15 S19-S12
      4 1833  837 (N09W59 L=255)  1427    245        3 (1800/ 770-1867/ 907) -33.7 W54-W67 N06-N12
      5  819  832 (S22E06 L=189)   265     69        0 ( 792/ 806- 853/ 872) -21.8 E09-E04 S23-S20
      6 1419  755 (S09W30 L=226)   653    158        2 (1402/ 673-1438/ 831) -53.1 W27-W33 S13-S05
      7 1001  757 (S21W06 L=202)   323    111        0 ( 947/ 739-1063/ 783)   7.1 W02-W10 S22-S20
      8 1362  702 (S14W29 L=225)   290     97        0 (1309/ 671-1390/ 727)   0.0 W25-W31 S14-S12
      9  596  636 (S41E17 L=179)   766    209        2 ( 499/ 618- 693/ 663)  10.0 E25-E08 S43-S37
     10 1437  453 (S24W45 L=241)  2469    868        3 (1122/ 391-1717/ 521)  16.1 W22-W74 S32-S15
     11 1223  254 (S43W46 L=242)   386     53        0 (1194/ 238-1250/ 271)  21.8 W43-W49 S44-S42
  

    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