Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  869 1867 (N43E49 L=184)   519     27        1 ( 850/1846- 888/1889) -45.0 E52-E44 N40-N45
      2 1230 1744 (N46E12 L=220)  1671    240        3 (1163/1667-1310/1794)  74.1 E15-E06 N37-N52
      3  433 1643 (N15E62 L=170)   502     81        1 ( 403/1626- 472/1662)  26.6 E65-E58 N13-N16
      4 1566 1411 (N33W23 L=255)   932    149        5 (1540/1341-1590/1486) -73.3 W20-W25 N28-N38
      5 1750 1342 (N36W41 L=274)  1643    299        2 (1674/1202-1807/1442) -33.7 W32-W48 N28-N41
      6 1322 1420 (N25W06 L=238)   215     58        1 (1295/1403-1359/1441)   0.0 W03-W08 N24-N26
      7  780  774 (S26E07 L=225)   407    154        0 ( 722/ 730- 846/ 809)  49.4 E09-E03 S30-S22
      8 1219  693 (S18W21 L=254)   635    276        0 (1133/ 658-1311/ 766)  20.0 W14-W26 S20-S14
      9  998  482 (S37W17 L=250)   449     75        0 ( 959/ 453-1044/ 514)  63.4 W15-W18 S39-S33
     10 1137  345 (S41W35 L=267)   327     49        1 (1111/ 328-1168/ 364)  76.0 W32-W36 S42-S38
     11 1320  229 (S41W61 L=293)   479     56        0 (1296/ 208-1347/ 255)  -5.7 W55-W65 S41-S39
  

    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