Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  401 1566 (N42E51 L=110)   700    146        3 ( 366/1541- 446/1588) -11.3 E55-E45 N40-N43
      2  663 1537 (N40E26 L=135)   567    133        4 ( 600/1509- 718/1572)  18.4 E30-E21 N37-N42
      3 1732 1417 (N23W57 L=218)   447     35        0 (1709/1397-1756/1440) -45.0 W54-W59 N21-N24
      4 1622 1275 (N16W43 L=204)   385    117        0 (1575/1218-1670/1325) -58.0 W39-W45 N11-N19
      5 1643 1191 (N10W43 L=204)   283     45        1 (1628/1167-1660/1218)  45.0 W41-W44 N09-N12
      6  922  913 (S02E07 L=154)   669    233        1 ( 901/ 802- 942/1037)  77.0 E08-E05 S08-N05
      7  312  664 (S16E55 L=106)   787    205        2 ( 261/ 633- 359/ 703)  22.6 E61-E49 S18-S13
      8 1032  606 (S22E02 L=159)   330     74        1 ( 997/ 590-1075/ 623)   0.0 E04-W01 S22-S20
  

    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