Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1100 1883 (N58E33 L=218)  1076     67        1 (1069/1846-1125/1925) -40.2 E44-E27 N53-N64
      2  740 1740 (N33E44 L=208)   412     80        0 ( 722/1710- 763/1780) -32.0 E47-E39 N31-N36
      3 1240 1679 (N44E05 L=247)   878    153        3 (1167/1630-1282/1719)  -5.2 E11-W00 N40-N47
      4  799 1633 (N27E32 L=219)   251     62        2 ( 779/1596- 815/1656) -63.4 E34-E30 N25-N29
      5 1609 1562 (N48W31 L=283)  1273    199        5 (1538/1479-1681/1614) -17.4 W23-W39 N42-N51
      6 1727 1432 (N42W44 L=296)   697     69        2 (1694/1399-1753/1463) -36.9 W39-W48 N39-N43
      7 1613 1254 (N25W32 L=284)  2080    510        5 (1466/1048-1749/1423) -29.7 W18-W47 N16-N32
      8  168 1007 (S24E65 L=186)  1647    200        3 ( 135/ 910- 206/1085)  70.0 E71-E59 S29-S18
      9 1079  623 (S25W14 L=266)  1313    285        2 ( 958/ 574-1218/ 664)  40.2 W08-W21 S31-S19
     10 1295  533 (S24W31 L=283)   308    135        0 (1239/ 491-1335/ 574)  76.0 W29-W32 S28-S20
     11  796  530 (S39E01 L=250)   233     70        0 ( 776/ 513- 826/ 559) -14.0 E03-W01 S40-S37
     12 1413  354 (S32W51 L=303)  2315    356        2 (1229/ 328-1547/ 385)  17.8 W36-W64 S37-S27
     13  968  164 (S63W42 L=294)  1945    154        2 ( 928/ 106-1024/ 232) -12.1 W28-W70 S67-S56
  

    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