Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  456 1668 (N59E55 L= 29)   715     83        0 ( 421/1640- 497/1691)   8.7 E62-E46 N56-N60
      2 1614 1416 (N19W49 L=133)   637    217        0 (1548/1355-1686/1451) -48.8 W43-W52 N13-N21
      3  634 1408 (N36E22 L= 62)   436    128        1 ( 574/1390- 695/1421) -16.7 E27-E17 N34-N37
      4  275 1060 (N17E54 L= 30)   526     56        2 ( 256/1028- 298/1093) -76.0 E55-E51 N15-N19
      5 1878 1006 (S11W64 L=148)   778     89        0 (1839/ 981-1911/1039)   5.2 W58-W69 S12-S09
      6 1045  568 (S22E06 L= 78)   967    270        1 ( 953/ 489-1129/ 669)  34.7 E13-W00 S26-S16
      7  567  570 (S16E38 L= 46)   478     91        1 ( 539/ 526- 607/ 605)  39.8 E41-E34 S18-S13
      8 1513  556 (S32W27 L=111)   371    113        0 (1471/ 511-1553/ 603)  35.5 W23-W30 S34-S29
      9 1565  506 (S37W33 L=117)   468     93        0 (1526/ 485-1609/ 542)  15.9 W28-W36 S38-S35
     10 1090  461 (S31E05 L= 79)   903    335        0 ( 947/ 396-1219/ 532)  14.0 E16-W04 S33-S27
     11 1350  445 (S37W14 L= 98)   267     72        0 (1318/ 421-1387/ 466)  11.3 W11-W16 S38-S35
     12  673  427 (S27E36 L= 48)   572     83        1 ( 646/ 389- 700/ 464) -78.7 E37-E33 S29-S24
  

    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