Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1465 1735 (N58W37 L=169)  1212    161        4 (1393/1710-1548/1754)  12.1 W25-W54 N53-N60
      2  921 1640 (N40E18 L=114)   386     79        2 ( 880/1607- 948/1661) -23.2 E21-E14 N37-N41
      3 1432 1623 (N46W27 L=159)   305     47        2 (1401/1605-1461/1635)   0.0 W23-W29 N45-N47
      4  664 1450 (N23E30 L=102)   488    154        0 ( 640/1391- 684/1519) -63.4 E32-E28 N19-N27
      5 1212 1464 (N31W07 L=139)   272     88        1 (1191/1422-1225/1509) -59.0 W05-W08 N28-N33
      6  488 1421 (N19E43 L= 89)   747    136        3 ( 446/1390- 552/1476) -15.3 E48-E37 N16-N21
      7 1744 1413 (N34W55 L=188)   960    141        1 (1682/1370-1795/1446) -13.0 W48-W61 N31-N35
      8  580 1413 (N19E35 L= 97)   306     81        2 ( 563/1366- 597/1445)  90.0 E36-E34 N16-N21
      9  462  816 (S19E36 L= 97)   796    151        4 ( 420/ 751- 516/ 858)  48.8 E39-E31 S23-S15
     10 1707  650 (S14W55 L=187)   479     35        2 (1688/ 621-1730/ 675) -53.1 W52-W57 S16-S12
     11 1410  412 (S34W40 L=172)   617    188        1 (1337/ 364-1492/ 453)  45.0 W36-W44 S38-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