Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1310 1763 (N53W39 L= 54)   854    103        3 (1268/1729-1360/1783)   0.0 W33-W46 N50-N55
      2  557 1482 (N36E34 L=341)   572    177        0 ( 515/1394- 598/1534)  73.3 E36-E31 N30-N40
      3  709 1434 (N32E21 L=354)   939    218        2 ( 644/1373- 792/1518) -34.7 E27-E14 N28-N38
      4 1539 1359 (N22W39 L= 54)   279    132        1 (1521/1304-1559/1414)  74.1 W36-W41 N18-N25
      5  859 1188 (N16E10 L=  5)   819    366        1 ( 721/1082- 987/1261) -34.5 E18-E01 N08-N20
      6 1258 1158 (N11W15 L= 30)   225     57        1 (1223/1144-1288/1174)   0.0 W13-W17 N10-N12
      7 1253 1088 (N07W15 L= 30)   393     63        1 (1213/1065-1289/1109)  14.0 W12-W16 N05-N08
      8 1597  817 (S12W37 L= 52)   448     84        1 (1560/ 786-1641/ 845)   9.5 W34-W40 S14-S10
      9  332  513 (S26E62 L=313)  1208    201        0 ( 272/ 433- 373/ 586) -84.8 E65-E59 S32-S21
  

    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