Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1196 1578 (N39W09 L=305)  7094   1374        3 ( 633/1406-1769/1746)  -8.4 E44-W59 N29-N47
      2  720 1196 (N10E20 L=275)  1774    427        3 ( 696/1001- 758/1399) -80.5 E23-E17 S01-N23
      3  419 1346 (N17E45 L=251)   753    141        3 ( 349/1301- 461/1397)  50.2 E50-E41 N14-N20
      4  177 1223 (N08E68 L=228)  2912    307        7 ( 154/1074- 203/1353) -55.5 E79-E62 S01-N15
      5 1557 1290 (N21W35 L=331)   486    128        0 (1518/1235-1577/1351)  58.0 W31-W36 N17-N25
      6 1478 1300 (N22W29 L=324)   282     73        0 (1459/1264-1503/1336) -63.4 W27-W30 N19-N23
      7  278 1228 (N09E55 L=240)   825    151        1 ( 258/1165- 306/1297) -76.0 E57-E51 N05-N13
      8  333  850 (S14E47 L=249)   937    241        2 ( 304/ 732- 354/ 961)  71.6 E51-E44 S21-S06
      9 1566  562 (S24W43 L=338)  3781    539        8 (1332/ 455-1775/ 670)  25.9 W26-W60 S32-S16
     10  708  632 (S25E19 L=277)   418    119        0 ( 656/ 600- 760/ 660)  35.5 E22-E15 S27-S22
     11 1189  415 (S37W17 L=313)   428     89        0 (1144/ 397-1232/ 434)  18.4 W14-W20 S38-S35
  

    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