Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  682 1803 (N55E32 L= 35)   641     43        0 ( 655/1777- 713/1823)  36.9 E35-E28 N52-N57
      2 1172 1506 (N27W12 L= 79)   266     72        0 (1151/1468-1201/1538)  51.3 W10-W14 N24-N29
      3 1426 1469 (N24W29 L= 96)   409     84        3 (1388/1448-1462/1499) -21.8 W26-W32 N22-N25
      4 1627 1436 (N21W45 L=112)   332     36        1 (1605/1413-1649/1461)  26.6 W42-W47 N19-N22
      5 1552 1386 (N18W37 L=104)   252     20        0 (1534/1368-1573/1408) -63.4 W35-W39 N17-N19
      6 1625 1365 (N17W43 L=110)   260     64        0 (1594/1341-1649/1391)  26.6 W39-W45 N15-N18
      7  261 1325 (N20E56 L= 10)   638    107        1 ( 239/1269- 291/1368)  80.5 E58-E53 N16-N22
      8 1587  872 (S13W36 L=103)   247     48        0 (1560/ 849-1607/ 889)  14.0 W34-W38 S14-S12
      9  388  558 (S28E52 L= 14)  1031    120        0 ( 339/ 527- 430/ 617)  19.7 E59-E45 S30-S24
     10  759  316 (S49E30 L= 37)   643     88        0 ( 729/ 289- 793/ 352) -90.0 E32-E25 S51-S45
  

    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