Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1107 1669 (N38E15 L=300)  1266    208        2 (1047/1578-1158/1756) -40.2 E21-E07 N32-N43
      2 1589 1663 (N57W30 L=344)  1077    117        1 (1536/1630-1644/1686)  27.6 W21-W44 N50-N63
      3  838 1553 (N22E26 L=288)  1128    354        4 ( 745/1398- 918/1679) -32.5 E38-E16 N14-N28
      4  285  939 (S27E50 L=265)  1502    338        1 ( 204/ 811- 393/1014)  36.9 E60-E39 S36-S20
      5  695  681 (S32E12 L=303)   781    255        1 ( 580/ 666- 825/ 697)  21.3 E20-E02 S35-S28
      6 1215  519 (S28W27 L=342)   737    151        1 (1184/ 451-1243/ 586) -60.9 W24-W29 S32-S23
      7 1484  476 (S22W48 L=  2)   297     72        1 (1455/ 451-1512/ 508)  90.0 W46-W49 S24-S19
      8 1531  446 (S22W54 L=  8)   627    141        1 (1498/ 397-1563/ 502) -74.1 W51-W56 S25-S18
      9  504  395 (S61E28 L=286)   853     72        2 ( 480/ 361- 528/ 441) -81.9 E33-E24 S63-S56
     10 1359  398 (S30W44 L=358)  1026    306        0 (1226/ 359-1471/ 435)  14.0 W34-W54 S34-S27
     11  618  327 (S63E11 L=303)   934    122        1 ( 580/ 291- 674/ 375) -18.4 E16-E02 S65-S58
     12  836  273 (S58W16 L=330)   377     69        0 ( 805/ 249- 866/ 295)   0.0 W12-W20 S59-S56
  

    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