Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1355 1788 (N54W49 L= 37)   501     59        1 (1332/1758-1378/1809)  21.8 W45-W55 N51-N56
      2  510 1522 (N40E40 L=308)   945    251        3 ( 404/1489- 617/1553)  -6.8 E53-E28 N37-N42
      3  848 1485 (N36E10 L=339)   335    117        0 ( 825/1425- 865/1536)  69.4 E11-E08 N31-N39
      4 1080 1463 (N32W07 L=356)  1087    273        3 ( 974/1416-1198/1530) -31.0 W00-W15 N28-N37
      5 1752 1364 (N19W59 L= 47)   381     70        1 (1735/1330-1773/1401)  68.2 W57-W61 N16-N21
      6 1131 1249 (N18W08 L=357)   630    274        1 (1061/1157-1226/1306)  33.7 W03-W15 N12-N21
      7 1649 1156 (N08W43 L= 32)   391     74        0 (1623/1117-1674/1192) -63.4 W41-W45 N05-N09
      8 1580 1099 (N05W37 L= 25)   464    108        0 (1516/1084-1625/1116)  -7.1 W32-W40 N03-N06
      9  475  491 (S27E47 L=302)   661    193        0 ( 454/ 406- 524/ 558)  38.2 E53-E39 S33-S22
     10  596  496 (S27E35 L=313)   598    181        1 ( 541/ 416- 639/ 553) -79.7 E38-E33 S33-S22
     11 1429  424 (S39W28 L= 16)   671    153        0 (1368/ 390-1491/ 470)  26.6 W22-W32 S41-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