Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1408 1821 (N47W52 L=360)   647     84        0 (1368/1794-1445/1847)  11.3 W44-W60 N45-N49
      2 1120 1620 (N35W14 L=322)   276     79        0 (1080/1604-1155/1633) -11.3 W11-W16 N33-N36
      3  452 1441 (N31E36 L=272)   762    170        0 ( 405/1377- 509/1510) -42.3 E42-E30 N26-N36
      4 1680 1375 (N13W48 L=356)   998    160        3 (1626/1317-1745/1431) -36.9 W44-W53 N09-N17
      5  717 1398 (N25E16 L=292)   208     80        0 ( 680/1376- 756/1421)  18.4 E18-E12 N23-N25
      6 1859 1011 (S10W59 L=  7)  2784    385       10 (1781/ 923-1924/1137) -72.9 W50-W67 S15-S02
      7  520 1052 (N06E31 L=277)  1753    404        6 ( 343/ 997- 700/1100)  -4.4 E44-E18 N02-N09
      8 1535  692 (S26W31 L=339)  4498    886        5 (1124/ 582-1901/ 830)   6.1 W01-W66 S30-S20
      9 1173  772 (S17W06 L=314)   180    111        0 (1132/ 742-1225/ 786)   9.5 W03-W09 S18-S16
     10  603  717 (S14E30 L=278)   510    230        1 ( 504/ 666- 684/ 739) -28.6 E36-E24 S18-S12
     11  934  630 (S23E10 L=298)   350    117        0 ( 904/ 599- 964/ 675)  51.3 E12-E07 S25-S20
     12  288  597 (S18E60 L=248)  1073    127        3 ( 261/ 539- 328/ 649)  54.5 E65-E54 S21-S14
     13 1792  525 (S41W65 L= 13)   611     41        2 (1765/ 507-1817/ 542)   7.1 W59-W71 S42-S39
     14 1642  477 (S42W46 L=354)   479     57        2 (1609/ 459-1676/ 501) -12.5 W41-W50 S43-S40
     15  421  412 (S31E58 L=250)  1636    205        3 ( 355/ 360- 518/ 440)  -9.5 E71-E48 S34-S28
     16  746  350 (S40E32 L=276)   413     94        0 ( 700/ 325- 794/ 374)   0.0 E37-E27 S41-S38
  

    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