Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  588 1669 (N61E17 L=248)   583    124        0 ( 552/1615- 631/1720)  37.9 E22-E10 N57-N65
      2 1168 1537 (N32W25 L=290)   424     77        0 (1133/1516-1218/1556) -31.0 W22-W28 N29-N33
      3  198 1130 (N31E62 L=203)   735    123        2 ( 150/1095- 254/1163) -17.4 E70-E54 N28-N33
      4 1587  768 (S23W28 L=293)   367     91        0 (1565/ 725-1612/ 799)  36.9 W26-W30 S25-S21
      5 1557  648 (S30W24 L=289)   553    104        0 (1512/ 618-1612/ 677)   0.0 W20-W28 S31-S28
      6  524  419 (S17E52 L=214)   686    133        2 ( 480/ 372- 564/ 476)  76.0 E53-E49 S20-S12
      7  686  382 (S23E41 L=224)   286    110        1 ( 663/ 350- 717/ 418)  68.2 E42-E38 S25-S20
      8  483  324 (S23E65 L=200)  3305    276        7 ( 388/ 254- 556/ 378) -51.3 E78-E54 S30-S19
  

    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