Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  774 1764 (N60E20 L=315)   352     76        0 ( 740/1747- 814/1775)   6.3 E24-E15 N58-N60
      2  641 1542 (N41E27 L=308)   259     41        1 ( 616/1525- 664/1557)  26.6 E29-E25 N40-N42
      3 1294 1471 (N31W23 L=359)   798    183        3 (1204/1445-1381/1499) -18.4 W16-W29 N28-N33
      4 1304 1265 (N18W20 L=355)   771    273        2 (1257/1150-1386/1349)  45.0 W15-W26 N10-N22
      5 1797 1214 (N09W59 L= 34)   574     76        0 (1768/1178-1816/1243)  53.1 W55-W61 N06-N10
      6 1742 1145 (N06W52 L= 27)   859    145        4 (1688/1120-1792/1194)   5.7 W46-W56 N04-N08
      7 1385  600 (S25W21 L=357)   720    265        0 (1274/ 541-1482/ 683)  25.0 W13-W28 S28-S20
      8  602  447 (S30E37 L=298)   837    195        0 ( 581/ 367- 637/ 541)  47.3 E43-E30 S36-S23
      9 1513  418 (S41W36 L= 11)   499     87        1 (1477/ 398-1547/ 440)   8.1 W32-W39 S42-S39
     10 1014  394 (S37E06 L=329)   267     71        0 ( 984/ 374-1059/ 416)  18.4 E08-E02 S38-S35
     11  560  304 (S42E52 L=283)   608    116        0 ( 518/ 267- 612/ 336) -63.4 E55-E49 S46-S39
  

    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