Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  791 1501 (N19E26 L=230)   441    166        1 ( 777/1425- 810/1569) -60.9 E28-E23 N14-N23
      2  686 1509 (N18E33 L=223)   247     83        0 ( 672/1473- 700/1545) -53.1 E35-E31 N15-N19
      3 1105 1458 (N24E06 L=251)   349     97        0 (1072/1413-1138/1492) -26.6 E08-E02 N22-N26
      4  664 1356 (N08E30 L=227)   434     89        1 ( 648/1308- 690/1398) -80.5 E31-E27 N05-N11
      5  227 1171 (S12E59 L=198)  1435    157        1 ( 190/1102- 267/1252)  78.7 E62-E56 S16-S06
      6 1351 1175 (N13W16 L=272)   228     70        0 (1312/1159-1390/1189)  36.9 W13-W18 N11-N14
      7  610 1026 (S12E25 L=232)   287     62        1 ( 584/ 995- 640/1057)  68.2 E26-E23 S14-S09
      8  557 1019 (S14E28 L=228)   300     67        2 ( 524/ 992- 582/1046)   0.0 E30-E26 S15-S12
      9  192  997 (S23E59 L=198)   595    120        1 ( 161/ 958- 208/1027) -26.6 E62-E56 S24-S20
     10 1218  744 (S14W18 L=275)   322     40        0 (1188/ 726-1247/ 763)  45.0 W16-W19 S15-S12
     11  575  720 (S31E22 L=235)   459     87        0 ( 548/ 697- 618/ 750)   0.0 E24-E18 S33-S29
     12 1429  711 (S11W32 L=289)   324    104        2 (1396/ 665-1476/ 747)  59.0 W31-W34 S14-S08
     13  348  688 (S39E42 L=215)   522     99        1 ( 323/ 651- 377/ 729)  56.3 E44-E38 S42-S36
     14  395  548 (S49E39 L=218)   710    150        0 ( 348/ 492- 453/ 607) -33.7 E43-E32 S51-S45
     15 1121  534 (S29W19 L=276)   380     82        0 (1072/ 517-1167/ 551)  18.4 W16-W22 S31-S27
     16 1639  465 (S19W60 L=316)   527     71        0 (1607/ 448-1666/ 485)  36.9 W57-W62 S20-S17
  

    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