Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1348 1494 (N29W32 L=276)   564    194        0 (1263/1471-1447/1511) -26.6 W25-W38 N25-N31
      2 1569 1411 (N19W45 L=289)   335    105        0 (1524/1390-1621/1430) -29.7 W41-W48 N16-N20
      3  294 1093 (N20E51 L=193)   594     75        0 ( 262/1067- 332/1122) -23.2 E54-E47 N18-N22
      4  428  870 (N05E41 L=203)   316     57        1 ( 406/ 842- 457/ 893)  21.8 E43-E38 N03-N06
      5 1502  570 (S31W25 L=269)  3313    581        4 (1282/ 439-1681/ 746)  22.1 W07-W40 S37-S22
      6  411  630 (S09E49 L=195)   720     83        2 ( 380/ 590- 443/ 670)  33.7 E51-E45 S11-S06
      7  401  450 (S22E60 L=184)   452     46        0 ( 378/ 434- 429/ 467) -26.6 E62-E57 S22-S20
      8  684  367 (S31E38 L=206)  1027    230        0 ( 642/ 280- 757/ 458)  34.4 E47-E28 S37-S24
      9 1059  429 (S32E09 L=235)   427     51        1 (1028/ 410-1089/ 452)  11.3 E11-E06 S33-S30
  

    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