Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1386 1522 (N33W33 L=356)   766    222        1 (1303/1482-1490/1570) -33.7 W27-W39 N29-N37
      2  787 1540 (N40E14 L=309)   296     58        1 ( 758/1518- 812/1563)  14.0 E16-E12 N38-N42
      3 1481 1341 (N20W34 L=357)   738    189        1 (1435/1254-1537/1407)  36.0 W29-W40 N15-N24
      4 1846 1191 (N06W64 L= 28)   698     57        2 (1822/1159-1870/1220)  20.6 W60-W68 N04-N07
      5 1437 1192 (N12W28 L=351)   224     67        0 (1416/1162-1465/1219) -63.4 W26-W29 N09-N13
      6 1661  856 (S11W42 L=  5)   683    196        1 (1629/ 765-1683/ 953)  80.5 W39-W44 S17-S05
      7  464  554 (S22E45 L=278)   476     65        1 ( 447/ 523- 487/ 588) -76.0 E46-E42 S23-S19
      8  940  489 (S29E10 L=313)   292     94        0 ( 928/ 441- 953/ 533)  81.9 E11-E08 S32-S25
      9  743  438 (S31E26 L=297)   712    193        0 ( 717/ 352- 788/ 511)  47.5 E31-E20 S38-S26
     10  674  294 (S43E40 L=283)   530     68        0 ( 645/ 266- 709/ 324) -68.2 E42-E37 S46-S40
  

    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