Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1032 1819 (N54W12 L=347)   499     75        1 ( 990/1798-1069/1840)   6.3 W07-W16 N51-N55
      2 1175 1806 (N50W25 L=360)   401     46        1 (1150/1779-1201/1828) -76.0 W23-W28 N47-N52
      3  611 1707 (N49E29 L=305)   544     47        1 ( 581/1683- 636/1736) -26.6 E32-E26 N46-N51
      4 1367 1568 (N30W30 L=  5)   380     87        3 (1344/1517-1397/1609)  50.2 W26-W33 N26-N32
      5  870 1496 (N29E06 L=329)   269     62        1 ( 834/1477- 902/1516)   0.0 E08-E03 N28-N30
      6  394 1447 (N31E43 L=291)   343     52        0 ( 371/1420- 421/1472)  21.8 E45-E40 N29-N32
      7  490 1408 (N27E34 L=301)   485     93        1 ( 447/1391- 549/1426)   0.0 E37-E28 N26-N28
      8  316  655 (S17E55 L=280)   366     91        1 ( 290/ 619- 336/ 693) -76.0 E56-E53 S18-S14
      9 1713  649 (S30W49 L= 24)   436     63        1 (1683/ 626-1745/ 677)  21.8 W46-W52 S31-S28
     10  423  613 (S20E47 L=288)   479     92        2 ( 399/ 573- 460/ 643)  33.7 E48-E42 S22-S18
     11  682  449 (S34E32 L=303)   491    117        1 ( 640/ 418- 738/ 478)  16.7 E36-E26 S35-S32
  

    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