Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  743 1577 (N29E22 L=224)   339     74        0 ( 723/1534- 766/1612)  90.0 E23-E21 N26-N32
      2  610 1485 (N23E30 L=216)   605    139        1 ( 543/1459- 675/1507) -10.3 E35-E24 N21-N24
      3  420 1353 (N14E42 L=205)   385     69        1 ( 401/1319- 439/1393)  90.0 E43-E40 N11-N16
      4  991 1282 (N11E03 L=243)   469    106        0 ( 948/1239-1044/1328) -39.8 E06-W00 N08-N13
      5 1851 1188 (N11W59 L=305)  3690    637        3 (1736/1039-1971/1313) -26.6 W47-W78 N03-N18
      6 1447 1251 (N11W25 L=272)   801    234        1 (1326/1209-1524/1312) -13.0 W17-W30 N09-N15
      7  293 1031 (S06E50 L=197)  1158    221        1 ( 272/ 920- 331/1144)  81.3 E51-E46 S12-N01
      8  725  920 (S12E18 L=229)   233    110        0 ( 667/ 898- 774/ 955)  23.2 E21-E14 S13-S09
      9  948  779 (S19E04 L=243)  1163    280        0 ( 851/ 670-1095/ 861)  38.2 E09-W05 S26-S13
     10 1757  641 (S22W56 L=303)  3779    601        5 (1588/ 532-1879/ 772)  37.7 W44-W68 S30-S13
     11  754  527 (S37E17 L=229)   562    148        1 ( 689/ 478- 813/ 577)  28.6 E23-E12 S40-S33
  

    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