Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1198 1494 (N25W06 L=319)   349    121        1 (1172/1436-1213/1551) -66.8 W04-W07 N21-N29
      2  455 1390 (N12E41 L=272)   978    184        2 ( 411/1331- 504/1464) -41.6 E46-E36 N08-N16
      3 1663 1187 (N12W40 L=352)  1857    441        6 (1546/1071-1807/1265) -18.4 W30-W52 N06-N15
      4 1130 1186 (N05W05 L=318)   250    103        1 (1101/1156-1179/1211)  26.6 W02-W07 N03-N06
      5 1239 1097 (N01W12 L=325)   182     72        1 (1213/1064-1262/1139)  90.0 W10-W13 N00-N03
      6  618  983 (S11E25 L=288)   244     75        1 ( 576/ 964- 655/1005)  21.8 E27-E22 S12-S09
      7  158  826 (S23E64 L=249)   523    149        2 ( 129/ 797- 179/ 852)  45.0 E67-E61 S24-S21
      8 1565  826 (S11W37 L=349)   336     79        1 (1530/ 800-1607/ 850)  -8.1 W33-W40 S11-S09
      9 1199  673 (S25W15 L=328)  1508    367        2 (1031/ 577-1391/ 802)  42.0 W05-W25 S32-S14
     10  355  630 (S35E47 L=266)   295     57        0 ( 331/ 606- 377/ 652) -45.0 E48-E44 S36-S33
     11  879  542 (S37E05 L=308)   583    126        3 ( 816/ 526- 945/ 568)  18.4 E09-W00 S38-S34
     12 1569  405 (S36W51 L=  4)   364     80        0 (1538/ 385-1597/ 430)   0.0 W47-W54 S37-S34
  

    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