Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  267 1482 (N10E70 L=206)  1263     89        3 ( 241/1437- 300/1520)  59.0 E78-E65 N07-N12
      2  300 1438 (N08E63 L=214)   538    108        1 ( 260/1392- 327/1480)  74.1 E64-E59 N04-N11
      3 1449 1376 (N28W20 L=296)   285     78        0 (1434/1332-1460/1410) -68.2 W18-W21 N25-N30
      4 1380 1260 (N19W17 L=293)   212     72        0 (1328/1245-1411/1273)  21.8 W13-W18 N18-N20
      5 1493 1031 (N08W28 L=305)   267     76        1 (1463/ 991-1513/1059) -21.8 W25-W30 N06-N09
      6  224  973 (S21E57 L=219)   484     65        0 ( 207/ 934- 243/1003)  78.7 E59-E55 S23-S18
      7 1771  974 (N11W52 L=328)   355     56        0 (1744/ 952-1792/ 999)  36.9 W49-W53 N09-N12
      8 1550  949 (N05W34 L=310)   356     59        1 (1524/ 919-1577/ 983)  76.0 W32-W35 N02-N06
      9  803  762 (S22E08 L=268)   207     42        0 ( 780/ 747- 832/ 783)  18.4 E09-E06 S22-S20
     10  543  629 (S36E25 L=251)   709    175        1 ( 474/ 595- 610/ 693)  -9.5 E31-E19 S38-S33
     11 1503  448 (S25W50 L=326)   986    152        2 (1440/ 403-1560/ 491)  60.3 W45-W54 S29-S21
  

    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