Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  792 1803 (N45E31 L=260)  1120    102        2 ( 744/1777- 855/1828)  -9.5 E37-E25 N42-N47
      2 1157 1699 (N40W02 L=293)   397    111        0 (1107/1681-1212/1718)  15.9 E01-W06 N38-N42
      3 1617 1619 (N42W43 L=335)   639    105        1 (1577/1594-1671/1644) -11.3 W38-W49 N39-N43
      4  459 1486 (N19E44 L=247)  1379    246        3 ( 391/1397- 551/1578) -28.8 E54-E34 N13-N24
      5  736 1519 (N23E24 L=267)   328    120        1 ( 716/1467- 759/1569) -54.5 E26-E21 N19-N26
      6 1346 1290 (N14W17 L=309)   301    158        0 (1300/1228-1382/1363) -58.0 W13-W19 N10-N18
      7 1722 1283 (N19W45 L=337)   808    196        1 (1648/1238-1808/1336) -12.1 W39-W54 N16-N21
      8 1649 1197 (N12W39 L=330)   516    135        0 (1607/1139-1689/1261) -40.6 W35-W42 N09-N15
      9 1419 1187 (N09W23 L=314)   205     70        1 (1388/1173-1463/1207)   0.0 W20-W25 N07-N09
     10 1874  871 (S03W64 L=356)   494     66        0 (1855/ 841-1893/ 907)  90.0 W62-W67 S05-S01
     11  664  709 (S27E21 L=270)  1312    418        0 ( 570/ 544- 788/ 874)  55.2 E30-E13 S39-S16
     12 1145  831 (S16W09 L=301)   195     49        0 (1117/ 814-1180/ 850)   0.0 W07-W11 S16-S14
     13 1773  591 (S21W65 L=356)  1186    143        4 (1717/ 557-1828/ 622)  24.4 W59-W74 S24-S19
  

    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