Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  972 1518 (N38W09 L=126)  1345    323        1 ( 937/1363-1015/1644)  48.6 W02-W17 N28-N45
      2 1195 1608 (N38W29 L=145)   368     57        0 (1171/1579-1227/1626)   0.0 W26-W31 N35-N38
      3 1420 1584 (N30W43 L=160)   675    127        1 (1372/1559-1488/1606)  -5.2 W38-W49 N28-N31
      4 1034 1289 (N22W06 L=123)   685    114        2 ( 984/1243-1088/1339) -56.3 W03-W08 N18-N25
      5  160 1128 (N27E68 L= 49)  1956    187        3 ( 127/1048- 194/1188)  48.8 E75-E62 N22-N31
      6  865  992 (N08E10 L=107)   451    198        3 ( 781/ 956- 946/1017)  -5.7 E15-E05 N07-N09
      7  263  908 (N13E55 L= 62)   680    196        2 ( 221/ 828- 290/ 964)  53.1 E59-E51 N08-N16
      8 1783  817 (S24W48 L=164)  2295    327        2 (1696/ 676-1843/ 960)  41.2 W38-W55 S30-S16
      9  214  857 (N11E62 L= 55)   411     92        0 ( 196/ 817- 240/ 901) -59.0 E63-E59 N08-N13
     10 1357  819 (S12W15 L=132)   361    120        2 (1321/ 785-1395/ 877)  26.6 W12-W18 S13-S09
     11  683  777 (S01E25 L= 91)   325    121        0 ( 670/ 720- 694/ 840)  63.4 E26-E23 S03-N03
     12 1420  723 (S19W17 L=134)  1155    412        4 (1306/ 649-1567/ 763)  -3.2 W09-W28 S23-S14
     13  873  676 (S10E16 L=100)   221     85        0 ( 861/ 644- 891/ 727)  80.5 E17-E15 S12-S06
     14 1040  628 (S16E08 L=109)   330    156        0 (1018/ 574-1077/ 688)  82.9 E09-E04 S20-S12
     15 1219  665 (S18W04 L=120)   323     56        0 (1182/ 642-1250/ 684) -36.9 W01-W05 S19-S16
     16  487  449 (S18E51 L= 66)  1338    231        0 ( 383/ 407- 559/ 516)   8.1 E62-E41 S20-S14
     17 1068  333 (S38E15 L=102)  4387    695        4 ( 796/ 207-1310/ 467)   9.7 E42-W05 S42-S33
     18 1564  352 (S54W30 L=146)   456     51        0 (1537/ 333-1599/ 370)   0.0 W25-W34 S54-S52
  

    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