Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  666 1573 (N50E10 L=191)   275     42        2 ( 643/1554- 695/1595) -11.3 E12-E07 N48-N51
      2 1551 1528 (N20W50 L=250)   414     91        3 (1509/1511-1597/1558) -14.0 W46-W54 N17-N21
      3 1519 1384 (N13W40 L=241)   393     51        0 (1493/1353-1548/1409)  11.3 W37-W42 N11-N14
      4 1573 1336 (N09W42 L=242)   326     95        1 (1543/1294-1606/1377)  38.7 W38-W44 N07-N11
      5  713 1277 (N29E13 L=188)   368    112        1 ( 667/1259- 774/1296) -29.7 E15-E08 N26-N31
      6 1011 1073 (N10W00 L=201)   318    104        0 ( 968/1036-1047/1104) -68.2 E01-W02 N07-N12
      7  742 1073 (N16E15 L=185)   269    140        0 ( 688/1038- 785/1099) -51.3 E18-E13 N13-N18
      8 1613  888 (S17W32 L=232)   318     73        1 (1576/ 859-1635/ 913) -36.9 W29-W33 S18-S15
      9 1226  710 (S16W04 L=204)   282     70        0 (1190/ 686-1253/ 737)  14.0 W01-W05 S17-S14
     10 1867  706 (S39W59 L=260)   630     62        2 (1838/ 682-1900/ 727)  -4.4 W54-W67 S40-S36
     11  264  620 (S01E65 L=136)  2290    286        5 ( 166/ 559- 375/ 676) -20.9 E79-E54 S05-N03
     12 1086  644 (S16E06 L=194)   242     85        1 (1054/ 623-1129/ 659)   0.0 E08-E03 S17-S14
     13 1061  569 (S20E10 L=191)   264     87        0 (1031/ 530-1080/ 603) -90.0 E10-E08 S22-S17
     14  434  375 (S19E63 L=137)  4723    523        9 ( 334/ 225- 568/ 547)  76.5 E73-E58 S32-S06
     15 1595  444 (S47W26 L=226)  1272    169        3 (1525/ 409-1667/ 494)   0.0 W18-W34 S49-S44
     16  732  322 (S28E41 L=159)  2550    391        8 ( 694/ 154- 773/ 488)  38.0 E61-E29 S42-S17
     17  461  422 (S16E56 L=144)   913    119        3 ( 432/ 375- 496/ 473)  81.9 E59-E52 S19-S12
  

    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