Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1474 1661 (N38W37 L=221)   369     44        0 (1454/1635-1497/1680)  31.0 W34-W39 N36-N39
      2  403 1583 (N35E48 L=137)   467     87        0 ( 364/1564- 440/1600)  18.4 E51-E44 N33-N36
      3  673 1576 (N33E24 L=160)   507    153        0 ( 599/1562- 748/1588)   0.0 E29-E18 N31-N33
      4 1580 1371 (N18W38 L=222)   851    284        1 (1466/1316-1707/1437) -19.4 W30-W47 N14-N22
      5  838 1357 (N17E11 L=173)   205     89        0 ( 796/1342- 869/1389)   0.0 E13-E09 N16-N19
      6  361 1309 (N16E45 L=139)   457    107        1 ( 327/1279- 401/1344) -23.2 E48-E41 N14-N18
      7  976 1156 (N05E03 L=182)   476    136        0 ( 910/1102-1020/1214) -49.4 E06-W00 N01-N08
      8  938  947 (S08E05 L=179)   568    179        1 ( 878/ 847- 995/1013)  55.0 E09-E01 S13-S03
      9 1547  932 (S09W33 L=217)   292     60        0 (1520/ 903-1574/ 965)  53.1 W31-W34 S10-S06
     10 1286  842 (S14W16 L=200)   632    360        2 (1176/ 781-1449/ 921)  16.4 W09-W26 S17-S09
     11  146  777 (S14E69 L=115)  1529    262        5 ( 120/ 706- 175/ 874)  78.7 E73-E65 S18-S08
     12  591  777 (S17E28 L=156)   284     64        1 ( 575/ 749- 616/ 812)  63.4 E29-E26 S18-S14
     13  323  782 (S15E49 L=135)   509     84        1 ( 294/ 755- 347/ 809)  31.0 E51-E46 S17-S13
     14 1809  770 (S18W57 L=241)   650     85        1 (1774/ 734-1840/ 802)  33.7 W53-W60 S20-S16
     15 1816  655 (S25W60 L=244)   614    154        1 (1770/ 638-1851/ 678)  14.0 W55-W63 S26-S23
     16  334  519 (S32E56 L=129)  1236    294        4 ( 314/ 424- 364/ 641)  60.3 E59-E49 S38-S24
     17 1199  517 (S35W12 L=196)   957    200        2 (1127/ 442-1278/ 599)  50.7 W06-W16 S40-S29
     18  615  438 (S39E34 L=150)  1531    197        4 ( 554/ 366- 664/ 527)  40.2 E41-E28 S44-S33
  

    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