Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  823 1729 (N38E36 L=142)   368     49        1 ( 799/1705- 851/1752)  -9.5 E39-E33 N36-N39
      2  863 1588 (N29E25 L=153)   706    167        1 ( 784/1551- 926/1626)  53.1 E29-E21 N24-N32
      3 1632 1271 (N26W36 L=214)  2393    749        1 (1421/ 997-1761/1544) -38.7 W15-W50 N12-N40
      4  545 1454 (N13E41 L=137)   629    103        1 ( 512/1405- 588/1498)  81.9 E42-E39 N09-N16
      5 1278 1362 (N24W08 L=186)   768    235        3 (1188/1300-1345/1448) -26.6 W00-W12 N20-N27
      6  545 1351 (N07E37 L=141)   256     49        0 ( 516/1327- 570/1367) -14.0 E39-E35 N06-N08
      7  238 1240 (S05E61 L=117)   485     79        1 ( 216/1203- 263/1273)  59.0 E63-E58 S07-S02
      8  455 1041 (S13E36 L=142)   712    180        4 ( 414/ 952- 525/1093) -83.7 E40-E31 S18-S09
      9  930  895 (S12E03 L=175)   323    106        1 ( 875/ 872- 980/ 907)   9.5 E06-W00 S12-S09
     10  660  823 (S22E18 L=160)   233     87        1 ( 624/ 803- 692/ 855)  36.9 E20-E16 S23-S18
     11 1662  723 (S05W49 L=227)   898    233        2 (1594/ 662-1720/ 765)  -8.7 W43-W56 S07-S02
     12 1484  690 (S11W37 L=215)   295     74        2 (1456/ 657-1514/ 719) -26.6 W34-W38 S13-S09
     13 1582  649 (S11W46 L=224)   659    134        4 (1544/ 603-1617/ 697) -80.5 W42-W47 S14-S08
     14 1053  629 (S25W11 L=189)   233     67        0 (1021/ 609-1084/ 653)  53.1 W09-W12 S26-S22
     15 1503  583 (S17W42 L=221)   955    219        2 (1402/ 560-1591/ 617)  30.3 W35-W48 S20-S13
     16  413  537 (S47E40 L=138)  1185    153        2 ( 351/ 499- 488/ 570)  24.2 E50-E30 S51-S42
     17 1566  509 (S20W52 L=230)   421     58        3 (1537/ 486-1592/ 534)  56.3 W48-W54 S21-S18
     18 1567  428 (S25W58 L=236)  1162    213        0 (1482/ 355-1658/ 487)  69.0 W54-W63 S31-S18
  

    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