Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  508 1496 (N09E44 L=280)  2896    552        2 ( 453/1260- 591/1706) -62.6 E50-E33 S03-N24
      2 1760 1472 (N42W44 L=  8)   511    102        2 (1738/1444-1788/1491)  35.5 W40-W47 N39-N44
      3 1642 1399 (N31W31 L=356)   329     50        0 (1608/1378-1666/1415) -11.3 W28-W33 N30-N32
      4 1295 1314 (N16W09 L=333)   455    130        0 (1262/1243-1323/1382) -49.4 W05-W11 N12-N19
      5  194 1191 (S13E62 L=262)  2500    502        4 ( 124/1146- 307/1269)  -7.4 E73-E50 S15-S08
      6 1342 1097 (N04W16 L=341)   205     51        1 (1312/1076-1375/1118)  14.0 W14-W18 N03-N05
      7 1714  974 (N08W44 L=  8)   373     68        0 (1695/ 936-1729/1015) -36.9 W41-W46 N05-N09
      8  723  937 (S19E16 L=309)   341    115        1 ( 667/ 921- 785/ 963)  29.7 E19-E11 S20-S16
      9 1422  708 (S15W32 L=356)   479    112        0 (1373/ 668-1485/ 758)  63.4 W29-W34 S18-S10
     10 1701  603 (S12W57 L= 21)   426     73        0 (1680/ 569-1715/ 632) -36.9 W53-W58 S13-S10
     11 1565  440 (S25W54 L= 19)  1662    202        4 (1505/ 407-1629/ 486)  15.3 W48-W59 S27-S21
  

    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