Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  440 1599 (N17E57 L=253)  2694    433        8 ( 338/1501- 512/1738)  90.0 E63-E53 N08-N27
      2  482 1481 (N10E47 L=263)   372    101        2 ( 457/1441- 503/1518) -59.0 E49-E44 N08-N13
      3  384 1469 (N08E55 L=255)   670    104        1 ( 351/1438- 426/1492)  15.9 E58-E51 N05-N09
      4  348 1353 (N00E53 L=257)   407     66        2 ( 318/1326- 374/1376)  45.0 E55-E51 S01-N02
      5  458 1080 (S13E36 L=274)   619    129        1 ( 425/1019- 494/1133)  76.0 E37-E33 S17-S09
      6  439  824 (S29E34 L=277)   296     65        0 ( 423/ 786- 463/ 846)  63.4 E34-E31 S31-S27
      7  977  752 (S20W04 L=314)   392    111        0 ( 926/ 722-1042/ 778)  39.8 W01-W07 S23-S17
      8  315  678 (S42E45 L=265)   946    264        1 ( 255/ 631- 381/ 718)  32.7 E52-E37 S46-S37
      9 1726  571 (S10W62 L= 12)   886    197        4 (1677/ 546-1774/ 598)  15.3 W56-W68 S11-S08
     10 1415  438 (S27W43 L=354)   583    169        0 (1342/ 415-1481/ 471)   9.5 W37-W49 S27-S25
     11 1365  264 (S39W52 L=  2)  1097    239        1 (1338/ 208-1385/ 329) -30.3 W47-W59 S42-S34
  

    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