Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  691 1520 (N18E30 L=246)   499     61        1 ( 652/1495- 735/1539)  18.4 E32-E26 N16-N19
      2 1715 1442 (N33W43 L=318)  1307    223        1 (1646/1344-1775/1515) -42.0 W36-W48 N27-N37
      3  370 1366 (N05E49 L=227)   624    100        2 ( 345/1310- 390/1414) -50.2 E51-E46 N01-N07
      4  154 1155 (S09E66 L=209)  1515    184        4 ( 118/1065- 176/1249)  90.0 E71-E63 S14-S03
      5 1847 1155 (N17W55 L=330)   897    153        2 (1778/1125-1915/1179)  13.2 W47-W64 N15-N20
      6  448 1080 (S11E37 L=238)   230     74        0 ( 424/1051- 466/1114)  63.4 E38-E35 S13-S09
      7 1895  981 (N09W63 L=338)  5252    258        6 (1840/ 869-1961/1103) -69.0 W54-W73 N01-N14
      8 1021  711 (S25W05 L=281)   459    141        0 ( 973/ 648-1070/ 767)  66.0 W03-W07 S29-S20
      9 1175  631 (S27W17 L=292)  1438    341        1 (1052/ 548-1337/ 680)  29.4 W10-W27 S33-S22
     10 1677  463 (S25W63 L=338)  4417    433        7 (1608/ 370-1766/ 610) -49.1 W54-W82 S30-S15
     11 1145  489 (S36W19 L=295)  1316    319        1 (1012/ 455-1293/ 537)  11.3 W09-W29 S38-S31
     12 1116  363 (S45W23 L=298)  3255    735        3 ( 821/ 252-1447/ 468)   1.8 E03-W60 S47-S41
  

    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