Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  450 1637 (N35E59 L= 66)   681    109        0 ( 417/1605- 502/1680)  66.8 E61-E55 N32-N39
      2 1122 1590 (N39W01 L=126)   954    176        2 (1053/1551-1214/1630) -18.4 E04-W08 N36-N41
      3  635 1490 (N27E33 L= 92)  1227    190        3 ( 570/1417- 697/1555) -26.6 E40-E26 N22-N31
      4  917 1377 (N23E10 L=114)   337    111        0 ( 878/1335- 950/1425)  56.3 E12-E08 N19-N25
      5 1227 1362 (N24W10 L=135)   336     88        3 (1208/1324-1249/1412) -78.7 W09-W12 N21-N27
      6  190 1031 (S06E63 L= 62)  2651    245        8 ( 138/ 924- 243/1118)  47.3 E71-E55 S13-N00
      7 1413 1007 (N04W24 L=149)   226     69        0 (1379/ 977-1443/1030) -26.6 W22-W26 N02-N04
      8  872  771 (S15E07 L=117)  2286    509        2 ( 696/ 620-1023/ 951)  53.1 E18-W00 S26-S02
      9 1692  743 (S10W49 L=174)   457     68        1 (1667/ 709-1708/ 784)  90.0 W47-W51 S12-S07
     10 1498  745 (S12W33 L=158)   380     76        2 (1478/ 701-1520/ 779)  78.7 W32-W35 S14-S09
     11 1803  688 (S13W64 L=188)   944     93        3 (1759/ 656-1838/ 715)  -4.1 W57-W71 S15-S11
     12 1474  599 (S21W35 L=159)   589    176        1 (1391/ 574-1529/ 633)  24.0 W29-W39 S23-S18
     13  616  587 (S30E25 L= 99)   676     94        2 ( 564/ 565- 660/ 618)   7.1 E29-E21 S31-S27
     14 1628  453 (S30W57 L=181)  2404    389        3 (1491/ 388-1757/ 569)  34.4 W44-W66 S35-S21
     15  369  483 (S41E54 L= 70)  1833    174        1 ( 315/ 422- 437/ 537)  18.9 E77-E42 S47-S35
     16  568  323 (S53E39 L= 85)   498     50        1 ( 544/ 307- 602/ 345)  15.3 E44-E33 S54-S50
  

    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