Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  625 1804 (N31E54 L=122)   790    137        2 ( 592/1784- 655/1835)   0.0 E57-E50 N29-N33
      2 1308 1661 (N40E01 L=176)   413     49        0 (1278/1630-1336/1688) -31.0 E03-W02 N37-N41
      3  655 1629 (N19E40 L=136)   486     91        0 ( 609/1606- 704/1647)  26.6 E43-E37 N17-N21
      4  747 1522 (N14E30 L=146)   578    181        1 ( 719/1452- 782/1606) -52.1 E34-E25 N10-N19
      5 1517  985 (N04W29 L=205)   398    111        1 (1465/ 966-1574/1002)  20.6 W25-W33 N01-N05
      6 1609  859 (N00W39 L=215)   380     79        1 (1572/ 834-1633/ 893)  56.3 W36-W40 S01-N02
      7  748  818 (S26E11 L=165)   629    183        1 ( 649/ 786- 831/ 854)  32.5 E17-E06 S29-S21
      8 1608  724 (S07W44 L=220)  1148    253        4 (1557/ 640-1648/ 812) -79.7 W40-W47 S12-S01
      9 1410  676 (S15W32 L=209)  1044    246        2 (1290/ 621-1519/ 749)  55.3 W27-W37 S21-S08
     10 1274  604 (S23W27 L=203)   698    158        2 (1204/ 579-1343/ 641)   5.7 W21-W31 S25-S20
     11 1038  440 (S40W20 L=197)  1811    367        5 ( 877/ 409-1213/ 474)  15.6 W08-W33 S44-S35
     12 1576  378 (S25W60 L=236)  1799    325        1 (1441/ 349-1658/ 412)  25.2 W49-W66 S29-S20
  

    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