Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1477 1735 (N58W39 L=172)  1438    210        3 (1388/1705-1580/1753)   5.2 W25-W60 N53-N60
      2  906 1642 (N40E19 L=114)   460     98        4 ( 856/1606- 938/1664) -23.2 E23-E16 N37-N41
      3 1420 1624 (N46W25 L=158)   332     44        1 (1391/1607-1444/1640)  18.4 W22-W28 N45-N47
      4  655 1452 (N23E31 L=102)   587    141        3 ( 623/1396- 692/1519) -63.4 E33-E28 N19-N27
      5 1201 1462 (N31W06 L=140)   265     75        1 (1183/1425-1215/1505) -68.2 W05-W07 N28-N33
      6  479 1426 (N19E44 L= 89)   848    146        3 ( 433/1393- 538/1467) -11.3 E48-E38 N16-N21
      7  568 1414 (N19E36 L= 97)   303     77        1 ( 553/1371- 583/1448) -68.2 E37-E34 N16-N21
      8 1728 1420 (N35W53 L=187)   867    135        1 (1665/1393-1784/1448)  -8.1 W46-W60 N33-N36
      9  232  878 (S18E58 L= 75)   393     75        2 ( 216/ 839- 255/ 909)  51.3 E60-E54 S20-S15
     10  443  802 (S20E37 L= 96)  1060    222        3 ( 400/ 693- 505/ 858)  50.7 E41-E32 S27-S15
     11 1700  649 (S14W55 L=188)   439     45        0 (1680/ 622-1723/ 669)  33.7 W52-W56 S16-S13
     12 1462  435 (S32W43 L=176)  1063    373        0 (1323/ 363-1615/ 505)  36.4 W34-W54 S38-S24
  

    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